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ABSTRACT 


For decades, analog magnetic tape recording was the most pop- 
ular method for recording music, but has been replaced over the 
past 30 years first by DAT tape, then by DAWs and audio inter- 
faces [1]. Despite being replaced by higher quality technology, 
many have sought to recreate a "tape" sound through digital ef- 
fects, despite the distortion, tape “hiss”, and other oddities analog 
tape produced. The following paper describes the general process 
of creating a physical model of an analog tape machine starting 
from basic physical principles, then discusses in-depth a real-time 
implementation of a physical model of a Sony TC-260 tape ma- 
chine. 


"Whatever you now find weird, ugly, uncomfortable, and nasty 
about a new medium will surely become its signature. CD distor- 
tion, the jitteriness of digital video, the crap sound of 8-bit - all 
of these will be cherished and emulated as soon as they can be 
avoided." -Brian Eno [2]. 


1. INTRODUCTION 


While analog magnetic tape recording (see fig.|1) is rarely used 
in modern recording studios, the sound of analog tape is still of- 
ten sought after by mixing and mastering engineers. To that end. 
several prominent audio plugin manufacturers including Wave 

Universal Audi?| and U-Hg"| have created tape emulating plug- 
ins. Unfortunately, the existing literature on analog tape emulation 
is somewhat lacking. While Arnadottir et at. and Valimaki et 
al. [4] describe the emulation of tape echo/delay devices, and Vali- 
maki et al [5] describe the emulation of disk-based audio recording 
media, we were unable to locate any existing research directly dis- 
cussing digital emulation of the magnetisation process, a gap in 
research that this publication intends to fill. That said, Kadis 
and Camras [6] discuss musical use of analog tape recorders in a 
useful technical manner, and Bertram [7] gives a in-depth techni- 
cal description of the physical underpinnings of analog magnetic 
recording; this work intends to build on their foundations. While 
tape machines also contain electronic circuits that contribute to 
the machine’s characteristic sound, this publication only consid- 
ers processes that relate directly to tape magnetisation. For read- 


https://www.waves.com/plugins/4j37-tape 
“https://www.uaudio.com/uad-plugins 


-in-bundles/magnetic-tape-bundle.html 
“https://u-he.com/products/satin 
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ers wishing to emulate tape machine circuits, a good overview of 
circuit modelling techniques can be found in [8]. 


Figure 1: An analog tape recorder (Sony TC 260) 


2. CONTINUOUS TIME SYSTEM 


Audio recorded to and played back from a tape machine can be 
thought of as going through three distinct processors: the record 
head, tape magnetisation, and the play head. 


2.1. The Record Head 


For an instantaneous input current J(t), the magnetic field output 
of the record head is given as a function of distance along the tape 
(x), and depth into the tape (y). Using the Karlqvist medium field 
approximation, we find [7]: 


H, (x,y) = = Ho(tan™ (@Al+=) 


+tan | (@P== *\) (1) 


(aia +¥) gy 


(9/2) + x)? + y? 


where H, and H, are components of the magnetic field H, g is 
the head gap, and Hp is the deep gap field, given by: 


1 
H,(#,y) = 5-Holn ( 


Hy = a (3) 


DAFX-1 


Proceedings of the 22™ International Conference on Digital Audio Effects (DAFx-19), Birmingham, UK, September 2-6, 2019 


where NV is the number of turns of wire around the head, and F is 
the head efficiency which can be calculated by: 


E= — (4) 


lAg dl 
1+ Urg Jcore A(l) 


where A, is the gap area, ju, is the core permeability relative to 
free space (410), and A(l) is the cross-sectional area at a point / 
along its circumferential length. 


2.2. Tape Magnetisation 


The magnetisation recorded to tape from a magnetic field can be 
described using a hysteresis loop, as follows [7]: 


M(a,y) = Froop(H (x, y)) (5) 


where F', oop is a generalized hysteresis function. 


Using the Jiles-Atherton magnetisation model, the following dif- 
ferential equation describes magnetisation in some direction (/) 
as a function of the magnetic field in that direction (H) [9]: 


dM __(1~0)6u(Man=M)_ | dMan 
dH  (1—c)dsk—a(Man—M) | ° dH 


where c is the ratio of normal and anhysteric initial susceptibilities, 
k is a measure of the width of the hysteresis loop, a is a mean field 
parameter, representing inter-domain coupling, and 65 and 64 are 
given by: 


1 if H is increasin 
bs = ee ibe (7) 
—1 if His decreasing 
1 if dg and Ma, — M have the same sign 
dm = (8) 
0 otherwise 
Man is the anisotropic magnetisation given by: 
A M 
Man = M,L(—+") 9) 
a 


where //, is the magnetisation saturation, a characterizes the shape 
of the anhysteric magnetisation and L is the Langevin function: 


Ta) nihG= - (10) 


2.3. Play Head 
2.3.1. Ideal Playback Voltage 


The ideal playback voltage as a function of tape magnetisation at 
a point x along the tape is given by [7]: 


6/2 M bog 
dy' h(x! +2,y')- M(a',y') 

5/2 dx 
(1) 
where JN is the number of turns of wire, W is the width of the 
playhead, E is the playhead efficiency, v is the tape speed, 6 is 
the thickness of the tape, and uo is the permeability of free space. 


Note that V(2) = V(vt) for constant v. h(a, y) is defined as: 


V(x) = NWEvpo de! / 


H(x,y) 
NIE 


where A(a, y) can be calculated by eqs. (1) and 2}. 


A(z,y) = (12) 


2.3.2. Loss Effects 


There are several frequency-dependent loss effects associated with 
playback, described as follows [1]: 


kd, [1 — a) [senta/?)) 
i= t 1 
V(t) = Volt)[e™]| = ae (13) 
for sinusoidal input Vo(t), where k is the wave number, d is the 
distance between the tape and the playhead, g is the gap width of 
the play head, and again 6 is the thickness of the tape. The wave 
number is given by: 


k= aut (14) 
v 
where f is the frequency and v is the tape speed. 
3. DIGITIZING THE SYSTEM 
3.1. Record Head 
For simplicity, let us assume, 
A(x,y,t) = H(0, 0,1) (15) 
In this case Hy = 0, and Hz = Ho. Thus, 
NEI(t 
iC eee (16) 
g 
or, ; 
H(n) = NEI(n) (17) 
g 


3.2. Hysteresis 


Beginning from eq. (6). we can find the derivative of / w.r.t. time, 


as in [9]: 


(1-0)5y(MsL(Q)-M)_ fF 
dM _ (a)Sgk-a(M.L(Q)—My 


dt 1—ca*®L'(Q) 
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(18) 
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where Q = , and 


L'(z) = 5 —coth?(r) +1 (19) 


Note that eq. can also be written in the general form for non- 
linear Ordinary Differential Equations: 


dM 


= ft M, a) (20) 


H 
where ti = [3 ; 


Using the trapezoidal rule for derivative approximation, we find: 


A(n) — H(n—1) 


7 H(n—1) (21) 


H(n) =2 
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We can use the Runge-Kutta 4th-order method to find an ex- 
plicit solution for M(n): 


kn =Tf(n 5 Mn 1) A ai(n *)) 
ks=TF(n 5 M(n 1) © ai(n 5)) (22) 
ka = T§(n,M(n he ka, t(n)) 

M(n) = M(n—1)4 2 = = = 


We use linear interpolation to find the half-sample values used 
to calculate kz and k3. Note that many audio DSP systems pre- 
fer lower-order implicit methods such as the trapezoidal rule to 
solve differential equations rather than a higher-order method like 
the Runge-Kutta method [8]. However in this case, it was found 
that the lower-order methods quickly became unstable for high- 
frequency input, particularly when the input is modulated by a bias 


signal (see Section|4.3). 


3.2.1. Numerical Considerations 


To account for rounding errors in the Langevin function for values 
close to zero, we use the following approximation about zero, as 


in [9]: 


iGa= coth(x) — +4 for |x| > 107+ (23) 
a 3 otherwise 
1 —coth?(x) +1. for|z| > 107+ 
3 otherwise 


Additionally, tanh(«), and by extension coth(x) is a rather com- 
putationally expensive operation. With this in mind, for real-time 
implementation, we approximate coth() as the reciprocal of a 
Gaussian continued fraction for tanh(«) , hamely 


tanh(x) = (25) 
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Figure 2: Digitized Hysteresis Loop Simulation 


3.2.2. Simulation 


The digitized hysteresis loop was implemented and tested offline 
in Python, using the constants M,, a, a, k, and c from (Uj. 
For a sinusoidal input signal with frequency 2kHz, and varying 
amplitude from 800 - 2000 Amperes per meter, fig. [2] shows the 
Magnetisation output. 


3.3. Play Head 
By combining eq. with eqs. and (16), we get: 
V(t) = NW EvpogM (t) (26) 


or, 


V(n) = NW EvpogM(n) (27) 


3.3.1. Loss Effects 


In the real-time system, we model the playhead loss effects with 
an FIR filter, derived by taking the inverse DFT of the loss effects 
described in eq. (13). It is worth noting that as in eq. (14), the 
loss effects, and therefore the FIR filter are dependent on the tape 
speed. 


The loss effects filter was implemented and tested offline in Python 
with tape-head spacing of 20 microns, head gap width of 5 mi- 
crons, tape thickness of 35 microns, and tape speed of 15 ips. The 
following plot shows the results of the simulation, with a filter or- 
der of 100. 


Tape Loss Effects vs. Frequency 


Amplitude [dB] 


103 104 
Frequency [Hz] 


Figure 3: Frequency Response of Playhead Loss Effects 


4. TAPE AND TAPE MACHINE PARAMETERS 


In the following sections, we describe the implementation of a 
real-time model of a Sony TC-260 tape machine, while attempt- 
ing to preserve generality so that the process can be repeated for 
any similar reel-to-reel tape machine. 


4.1. Tape Parameters 


A typical reel-to-reel tape machine such as the Sony TC-260 uses 
Ferric Oxide (y F203) magnetic tape. The following properties of 
the tape are necessary for the tape hysteresis process eq. (18): 
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e Magnetic Saturation (/,): For Ferric Oxide tape the mag- 
netic saturation is 3.5e5 (A/m) [12]. 


e Hysteresis Loop Width (4): For soft materials, k can be 
approximated as the coercivity, H, . For Ferric Oxide, 
H is approximately 27 kA/m [12]. 


e Anhysteric magnetisataion (a): Knowing the coercivity and 
remnance magnetism of Ferric Oxide (12, we can calculate 
a = 22 kA/m by the method described in : 


e Ratio of normal and hysteris initial susceptibilities (c): From 
,c=1.7e-1. 


e Mean field parameter (a): From , a= 1.6e-3. 


4.2. Tape Machine Parameters 
4.2.1. Record Head 


To determine the magnetic field output of the record head using 
eq. (7), the following parameters are necessary: 


e Input Current (I(n)): For the Sony TC-260 the input cur- 
rent to the record head is approximately 0.1 mA peak-to- 
peak [14]. 


e Gap Width (g): The gap width for recording heads can 
range from 2.5 to 12 microns [I]. 


e Turns of wire (V): The number of turns of wire is typically 
on the order of 100 [7]. 


e Head Efficiency (£): The head efficiency is typically on the 
order of 0.1 [7]. 


These values result in a peak-to-peak magnetic field of approxi- 
mately 5e5 A/m. 


4.2.2. Play Head 


Similar to the record head, the following parameters are needed 
to calculate the output voltage using eqs. and (note that 
values are only included here if notably different from the record 
head): 


e Gap Width (g): The play head gap width ranges from 1.5 to 
6 microns[]]. 


e Head Width (W): For the Sony TC-260, the play head 
width is 0.125 inches (note that this is the same as the width 
of one track on the quarter-inch tape used by the machine) 


; 


e Tape Speed (v): The Sony TC-260 can run at 3.75 inches 
per second (ips), or 7.5 ips [14]. Note that many tape ma- 
chines can run at 15 or 32 ips [1]. 


e Tape Thickness (5): Typical tape that would be used with 
the TC-260 is on the order of 35 microns thick ‘ 


e Spacing (d): The spacing between the tape and the play 
head is highly variable between tape machines. For a typi- 
cal tape machine spacing can be as high as 20 microns [I]. 


2 kHz Sine Wave Biased at 50 kHz 
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Figure 4: Example of a biased signal 


4.3. Tape Bias 


A typical analog recorder adds a high-frequency "bias" current to 
the signal to avoid the "deadzone" effect when the input signal 
crosses zero, as well as to linearize the output. The input current 
to the record head can be given by [6]: 


Theaa(n) = Tin (n) + Bcos(2r foiasnT) (28) 


Where the amplitude of the bias current B is usually about one or- 
der of magnitude larger than the input, and the bias frequency foias 
is well above the audible range. Figure/4]shows a unit-amplitude, 2 
kHz sine wave biased by a 50 kHz sine wave with amplitude 5. To 
recover the correct output signal, tape machines use a lowpass fil- 
ter, with a cutoff frequency well below the bias frequency, though 
still above the audible range [I]. 


For the Sony TC-260, the bias frequency is 55 kHz, with a gain 
of 5 relative to the input signal. The lowpass filter used to recover 
the audible signal has a cutoff at 24 kHz, though note that due to 
the frequency response of the playhead loss effects, the effects of 
this filter may be essentially neglible to the real time system. 


4.4. Wow and Flutter 


Each tape machine has characteristic timing imperfections known 
as “wow” and/or “flutter.” These imperfections are caused by mi- 
nor changes in speed from the motors driving the tape reels, and 
can cause fluctuations in the pitch of the output signal. To charac- 
terize these timing imperfections, we use a method similar to BB}: 
We recorded a pulse train of 1000 pulses through a TC-260, then 
recorded the pulses back from the tape. Figure[5]shows a section 
of a superimposed plot of the original pulse train against the pulse 
train recorded from the tape machine. From this data, we were able 
to generate a periodic function that accurately models the timing 
imperfections of the TC-260. The process was performed at both 
7.5 ips and 3.75 ips. In the real-time system, the timing imperfec- 
tion model is used to inform a modulating delay line, to achieve 
the signature "wow" effect of an analog tape machine. 
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Figure 5: Input pulse train superimposed with pulse train recorded 
from TC-260 


Record Head Constants 


Turns of wire (N) 100 
Head Efficiency (E) 0.1 
Record Head Gap (g) 6 (microns) 
Tape Constants 
Magnetic Saturation (V/;) 3.5e5 (A/m) 
Hysteresis Loop Width (k) 27 (kA/m) 
Anhysteric Magnetisation (a) 22 (kA/m) 
Ratio of magnetic susceptibilities (c) 1.7e-1 
Play Head Constants 
Play Head Width (W) 0.125 (in) 


Figure 6: Constant values for model implementation 


5. REAL-TIME IMPLEMENTATION 


We implemented our physical model of the Sony TC-260 as a VST 
audio plugin using the JUCE framework. Figure[7|shows the signal 
flow of the system in detail. We allow the user to control param- 
eters in real-time including the tape speed, bias gain, gap width, 
tape thickness, tape spacing, and flutter depth. 


Bias modulation is implemeted using eq. (28), where bias gain B 
and bias frequency fics are controlled by the user. 


The record head transfer function calculates the record head mag- 
netic field H from the input current J, and is implemented using 
eq. (17), with constant values shown in fig. [6] 


The hysteresis process calculates the tape magnetisation / from 
the record head magnetic field H, and is implementated using the 
Runge-Kutta method described in eq. 22}, with constant value de- 
fined in fig. [6] 


The flutter process is implemented using a modulated delay line 
as described in section with user controlled modulation depth. 


The play head transfer calculates the play head voltage V from 
the tape magnetisation /, using eq. and the loss effects filter 


described in section The gap width g, tape speed v, tape 
thickness 6, and spacing d are controlled by the user, and other 
constant values are shown in fig. [6] 


C/C++ code for the full real-time implementation is open-source 
and is available on GitHub/’] 


5.1. Oversampling 


If no oversampling is used, the system will be unstable for input 
signal at the Nyquist frequency, due to limitations of the trape- 
zoid rule derivate approximation used in eq. (21). To avoid this 
instability, early versions of the real-time implementation used a 
lowpass filter with cutoff frequency just below Nyquist. However, 
due to aliasing caused by the nonlinearity of the tape hysteresis 
model, oversampling is necessary to mitigate aliasing artifacts [8]. 
Further, the system must be able to faithfully recreate not only the 
frequencies in the audible range but the bias frequency as well. 
Since the TC-260 uses a bias frequency of 55 kHz and the 
minimum standard audio sampling rate is 44.1 kHz, a minimum 
oversampling factor of 3x is required. However, since the biased 
signal is then fed into the hystersis model, even more oversam- 
pling is required to avoid aliasing. With these considerations in 
mind, our system uses an oversampling factor of 16x. 


5.2. Results 


In subjective testing, our physical model sounds quite convincing, 
with warm, tape-like distortion, and realistic sounding flutter. The 
high-frequency loss and low-frequency “head bump” change cor- 
rectly at different tape speeds, and are approximately within the 
frequency response specifications of the TC-260 service manual 
. When the input to the plugin is silent, the hysteresis process- 
ing of the bias signal produces a very accurate “tape hiss” sound. 
The distortion and frequency response characteristics of our model 
are subjectively very close when compared to the output of an ac- 
tual TC-260, though not nearly close enough to “fool” the listener. 
Additionally, as the bias gain is lowered, the “deadzone” effect 
appears exactly as expected [6]. The largest difference between 
the model and the physical machine is the subtle electrical and 
mechanical noises and dropouts present in the physical machine, 
presumably caused by the age and wear-and-tear of the machine, 
which we did not attempt to characterize in our model. Figure[8| 
shows the results of tests performed on the real-time system, in- 
cluding an example of the “deadzone” effect, and the timing irreg- 
ularities or “flutter”. Figure [9] shows a comparison of hysteresis 
characteristics between the real-time software model and a phys- 
ical Sony TC-260 tape machine. Note that some differences be- 
tween the two hysteresis loops may be due to the circuitry of the 
tape machine that we did not attempt to model in the real-time 
system. Audio examples from the real-time system can be found 
online] 


5.3. Evaluation 


While there is an audible difference between the real-time soft- 
ware model and a physical Sony TC-260, the most fundamental 
aspects of the tape machine sound including tape saturation, tape 


Shttps://github.com/jatinchowdhury18/ 
AnalogTapeModel 


https://ccrma.stanford.edu/~jatin/420/tape/ 
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Bias Signal 


Flutter (7) Tape Speed (v) 


Figure 7: Flowchart of realtime system: M is the oversampling factor, Hrecora is the transfer function of the record head, and Hptay is the 


play head transfer function including loss effects, and a de-biasing filter. 
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Figure 8: Testing results for real-time system: sine wave output 
with no biasing (above), input vs. output pulse train comparison 
(below). 


hiss, flutter/wow, and frequency response, are clearly audible and 
sound very accurate. The main distinctions between the two sys- 
tems can be attributed to the tape machine circuitry (in particular 
the TC-260 contains two shelving filters), as well as mechanical 
wear of the system, both elements that were not considered in our 
model. 


In our opinion, the strongest proof of the efficacy of our model 
is that the model responds accurately to the adjustement of model 
parameters. In particular, the hysteresis process reacted correctly 
to changes in input gain (saturating for overdriven input, or fading 
into tape hiss for underdriven input), as well as bias gain (satura- 
tion for overbiasing, or “deadzone” effect for underbiasing). Ad- 
ditionally, adjusting the loss effect parameters correctly demon- 
strated known tape machine phenomena including head “bump” 


Real-time System Hysteresis Loop 
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Figure 9: Test results comparing real-time system to Sony TC-260 
physical unit: hysteresis loop for real-time system (above), hys- 
teresis loop for TC-260 (below). 


(a resonance at the wavelength of the play head gap width), and 
spacing loss (filtering due to the spacing between the the play head 
and tape). The reader is invited to download the plugin (available 
with the source code) and evaluate the model for themselves. In 
conclusion, we believe that our model successfully approximates 
the physical tape recording process, however for those wishing to 
model a full tape machine, we suggest using this model in combi- 
nation with a model of the tape machine’s circuits. 


6. FUTURE IMPROVEMENTS 


6.1. Spatial Magnetic Effects 


The most obvious improvement to be made for the physical model 
is the inclusion of spatial effects of the tape. In particular, the ap- 
proximations made in eq. (U5), negate any effects caused by mag- 
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netisation along the longitudinal length of the tape, and into the 
depth of the tape. Including spatial effects would involve deriv- 
ing digital analogues for eqs. (Ip, and (Ut), and re-deriving 
eq. (22) to take an 2-dimensional magnetic field input at every 
timestep, rather than the zero-dimensional input it currently takes. 
This change would greatly increase the computational complexity 
of the system. At an oversampling rate of 16x, using just 100 spa- 
tial samples would be 1600x more computationally complex than 
the current system. 
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